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ABSTRACT 

Context. Planetary nebulae in the nearest large elliptical galaxy provide light element abundances difficult or impossible to measure 
by other means in a stellar system very different from the galaxies in the Local Group. 

Aims. The light element abundance pattern from many planetary nebulae (PNe) at a range of radial distances was measured from 
optical spectroscopy in the elliptical galaxy NGC 5128, which hosts the radio source Centaurus A. The PN abundances, in particular 
for oxygen, and the PN progenitor properties are related to the galaxy stellar properties. 

Methods. PNe in NGC 5128 covering the upper 4 mag. of the luminosity function were selected from a catalogue. VLT FORS1 
multi-slit spectra in blue and red ranges were obtained over three fields at 3, 9 and 15' projected radii (4, 8 and 17 kpc, for an adopted 
(-H . distance of 3.8 Mpc) and spectra were extracted for 51 PNe. 

Accurate electron temperature and density diagnostics are usually required for abundance determination, but were not available for 
most of the PNe. Cloudy photoionization models were run to match the spectra by a spherical, constant density nebula ionized by a 
black body central star. He, N, O and Ne abundances with respect to H were determined and, for brighter PN, S and Ar; central star 
luminosities and temperatures are also derived. 
Results. Emission line ratios for the 51 PNe are entirely typical of PN such as in the Milky Way. The temperature sensitive 
[O IIIJ4363A line was weakly detected in 10 PNe, both [O II] and [O III] lines were detected in 30 PNe, and only the bright 
[O III]5007A line was detected in 7 PN. For 40 PNe with Cloudy models, from the upper 2 mag. of the [O III] luminosity func- 
tion, the most reliably estimated element, oxygen, has a mean 12+Log(0/H) of 8.52 with a narrow distribution. No obvious radial 
gradient is apparent in O/H over a range 2-20 kpc. Comparison of the PN abundances with the stellar population, from the spectra 
f 1 ^ , of the integrated stellar light on the multi-slits and existing photometric studies, suggests an average metallicity of [Fe/H]=-0.4 and 

Q ■ [O/Fe]=0.25. 

Conclusions. The masses of the PN central stars in NGC 5128 deduced from model tracks imply an epoch of formation even more 
recent than found for the minority young population from colour magnitude studies. The PN may belong to the young tail of a recent, 
minor, star formation episode or derive from other evolutionary channels, perhaps involving binary stars. 

m . 

Key words. Galaxies: elliptical and lenticular, cD - galaxies: individual: NGC 5128 - galaxies: abundances - planetary nebulae: 
. general - stars: abundances 

1 . Introduction method (Jacoby et al. [T992l Ciardullo et al l2002l . The radial ve- 

locity, simply measured from one emission line, provides a kine- 

5^ , Planetary nebulae provide a multi-facetted probe of the galactic matic test particle for study of the galaxy pot e ntia i and in intra . 

environment. As a short-lived product of the evolution of low cluster gas of the cluster pote ntial and turbulence (Arnaboldi et 

mass stars, they sample the bulk, by both number and mass, of al [20031 . The abundances of a number of light elements can be 

the stellar population of galaxies of all Hubble types. Their size measured from emission line ratios in PN spectra, such as He 

makes them unresolved from the ground at distances greater than and N an ^ in particu i ari the a-elements O, Ne, S and Ar, which 

~ 1 Mpc but their strong emission lines makes them easy targets ^ difficult t0 study in all but high resolution and high signal-to- 

for observational work, even when observed against the high no j se spectra of individual stars 
stellar background of a galaxy bulge. These advantages have 

been exploited in a number of distinct areas. Statistics of PN As the closest example of a large early-type galaxy, 

in galaxies of various types show a late time indicator of the star NGC 5m (Hubble type SOp) occupies a central place in studies 

formation rate through the luminosity specific PN frequency <j> ( of reso i V ed stellar populations in a galaxy different to the Milky 

Ciardullo et al.[1989). The PN luminosity function (PNLF), de- Way and M31 ms large elliptical galaxy in the Centaurus 

spite theoretical difficulties, still proves to be a powerful distance group shows signs of major activity with an active nuc i eus , 

Faranoff-Riley Class I radio lobes, presence of dust and young 

Send offprint requests to: J. R. Walsh stars in its inner region (Graham 1979"i, a young (~0.3 Gyr) tidal 
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stream (Peng et al. 2002) and stellar shells in its outer regions 
(Malin et al. 1983); see Israel (119981 ) for an earlier review. This 
activity can be attributed to a minor merger which has had lit- 
tle influence on the bulk of the older passively evolving popula- 
tion (Woodley 2006), making NGC 5128 a nearby exemplar of 
more distant massive early-type galaxies; Harris (2010 1 provides 
a recent review of the underlying galaxy properties. The spec- 
troscopy of globular clusters by Beasley et al. (2008 ) indicates 
typical ages of 7-8 Gyr from stellar population fitting, reinforc- 
ing the presence of a large-scale intermediate age star formation 
episode. Fits to deep colour-magnitude diagrams also provide 
evidence for a minority, much younger, population of age 2-4 
Gvr (Reikuba et al. EOTTI ). 

As representative of the low mass stars, planetary nebulae in 
NGC 5128 have been catalogued over a period of two decades, 
beginning with the catalogue of Hui et al. d!993al ). 785 PN were 
discovered by [O III]5007A emission line and off-band imag- 
ing; from the PN [O III] magnitudes Hui et al. d!993bl ) deter- 
mined a PNLF distance of 3.5 Mpc. Independent measurements 
from the Mira period-luminosity relation and the luminosity of 
the tip of the red giant branch (Rejkuba et al. 2005 ), surface 
brightness fluctuations (Tonry et al. 12001), globular cluster lu- 
minosity function (Harris et al. 1988 ) and 42 classical Cepheid 
variables (Ferrarese et al. 12007) result in distance estimate in 
the range 3.4 to 4.1 Mpc. Harris et al. (2010) present a compre- 
hensive review of distance estimates to NGC 5128, and recom- 
mend a best-estimate of 3.8 Mpc, which is adopted here (dis- 
tance modulus 27.90mag.). Peng et al. (2004) extended the orig- 
inal catalogue of PN and found a further 356 PN by filter imag- 
ing; 780 out of a total of 1 141 were spectroscopically confirmed. 
Further emission line mapping and follow-up intermediate dis- 
persion spectroscopy has extended this list to over 1200 con- 
firmed PN (Rejkuba & Walsh 120061 ) . This number of PN makes 
NGC 5128 a rich source for statistical extra-galactic PN studies, 
rivalled only by the Milky Way and M 31 (Merrett et al. 120061 

Since there are both a large number of PN and they have been 
catalogued to large galactocentric distances (to 80 kpc along the 
major axis by Peng et al. (120041 )'). the radial velocities allow the 
dynamics of the halo mass distribution to be studied. Hui et al. 
( 119951 ) measured [O III] radial velocities for 431 out of their 
785 catalogued PN. The offset of the rotation axis from the mi- 
nor axis was attributed to evidence of triaxiality; fitting the ro- 
tation curve and velocity dispersion revealed a radially increas- 
ing mass-to-light ratio and hence the presence of a dark mat- 
ter halo. Peng et al. (2004) extended the PN kinematic work to 
780 PN and showed the large rotation along the major axis but 
with a pronounced zero-velocity twist produced by the triaxial- 
prolate mass distribution. The kinematics of the large popula- 
tion of globular clusters (with 563 available radial velocities - 
Woodley et al. BOTUl Woodley et al. 120071 Beasley et al. 120081 
show similar kinematics to the PN but with small differences, 
such as lower rotation amplitude (Woodley et al. 120071 ). 

The NGC 5128 globular clusters (GCs) also provide funda- 
mental evidence of the star formation history. Of the 605 con- 
firmed GCs in NGC 5128 (Woodley et al. l2U07l Woodley et al. 
120101 ). more than half have metallicity measurements and they 
divide roughly half and half into metal rich and metal poor above 
and below [Fe/H]= -1.0. The metal poor GCs ([Z/H]~ -1.3) 
are older with ages similar to Milky Way GCs and the metal rich 
GCs (Peng et al. 120041 Beasley et al. 2008]) have intermediate 
ages (4-8 Gyr) and metallicity ([Z/H]~ -0.5). There are many 
more GC candidates (e.g. Harris et al. 2004) and the total num- 
ber of GCs has been estimated at around 1300 (Harris 2010 1. 



NGC 5128 is close enough that individual stars can be re- 
solved with HST imaging (and from the ground with adap- 
tive optics in the near-IR) and several studies have been ob- 
tained outside of the bright bulge where crowding is lower. Soria 
et al. ( 119961 ) detected the Red Giant Branch (RGB) and the 
Asymptotic Giant Branch (AGB) from intermediate age stars. 
The same field was followed up with NICMOS photometry and 
the IR colour magnitude diagram shows stars well above the tip 
of the red giant branch for old stars, confirming the presence 
of intermediate age stars (Marleau et al. 12000 ). Colour magni- 
tude diagrams from HST WFPC2 were constructed for two halo 
fields at projected radii of 19 and 29 kpc (Harris et al. 1999 
Harris et al. |2p00). The colour magnitude diagrams are domi- 
nated by old stars. Under the assumption that the ages of the 
stars in NGC 5128 are similar to those of old globular clusters 
in our Galaxy, the mean metallicity of the two fields was found 

to be similar at [Fe/H] 0.4 with about one third of the stars 

in a metal poor component and the rest metal rich. The metallic- 
ities in the halo fields were compared by Harris et al. (2002) to 
the metallity distribution function (MDF) in a field at 7.4kpc, 
containing a mixture of outer bulge and inner halo. A broad 

MDF was found as in the halo fields, peaking at [Fe/H] 0.4; 

but subtracting the MDF of the halo fields reveals a peak at 
slightly higher metallicity (~ -0.2). Much deeper ACS imaging 
by Rejkuba et al. (120051) in a halo field at 36 kpc reached to the 
horizontal branch (core helium burning population) as well as 
showing the red giant branch, red clump and AGB bump. Again 
a broad MDF was found with an average metallicity of -0.64, 
but with a broad tail to higher metallicity (> 0). The age sen- 
sitive indicators imply an average age for the halo of 8^ Gyr. 
Comparison of the MDF in the four fields shows the peak shift- 
ing to lower metallicity with increasing projected radius but the 
distribution is similarly wide at all radii. 



In order to determine light element abundances of PN in 
NGC 5128 fairly high signal-to-noise spectroscopy is required. 
Close to the inner bulge the stellar continuum surface brightness 
is too high and only the strongest few lines can be measured. 
However in the outer regions by selecting bright PN, diagnos- 
tic lines of He, Ne, Ar and S as well as the brighter lines of O 
and N can be measured and abundances derived. Walsh et al. 
d 19991) performed deep ESO 3.6m long slit spectra of a few se- 
lected PN with a long slit spectrograph and measured lines in 
addition to the brightest line ([O III]5007A) in five PN. O/H 
could be determined in two PN and values of 12 + Log(0/H) 
~8.5 (i.e. [O/H] - 0.2, adopting the Solar oxygen abundance of 
8.69 from Scott et al. (f2009)) were derived. This initial work 
has been expanded to deeper spectra of many more PNe with 
the ESO VLT and FORS1 instrument in multi-slit mode. In sec- 
tion 2 the observations are described and the reduction of the 
spectra and derivation of the abundances are presented in sec- 
tion 3. The results for the PN line fluxes and derived abundances 
both for the individual PN and for summations of many spectra 
are presented in section 4. Since the diagnostic lines for elec- 
tron temperature and density, important for abundance determi- 
nation, are weak or undetected, photoionization modelling of 
all well-detected lines was undertaken using Cloudy and is de- 
scribed in section 5. Consideration of the lack of a gradient in 
the PN O abundances, comparison of the results with the stel- 
lar, photometrically-derived metallicity, and relation of the PN 
to the stellar populations in NGC 5128 are discussed in section 
6. Conclusions are collected in section 7. 
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2. Observations 

2.1. FORS imaging and spectroscopy 

The European Southern Observatory Very Large Telescope 
(VLT) FORS1 instrument mounted on VLT Unit Telescope 1 
was used for the observations described. FORS 1 has a number 
of modes for imaging; for spectroscopy, long slit, multi-slit, or 
multi-slit mask are available (see Appenzeller et al.|T998 for de- 
tails). The multi-object mode (MOS) allows 19 slitlets of height 
varying from 20 to 22 " to be placed over a 6.8x6.8' field of 
view. In the regions of NGC 5128 where PNe were catalogued 
by Hui et al (1993a) and ( 1993b ), the surface density is such that 
there are usually enough PNe to be able to fill the slitlets, thus 
ensuring an optimal match between target density and instru- 
ment multiplex. This match was realised in practice for regions 
over the bright disk of NGC 5128, whilst in the outer regions, 
around 50% of the slitlets could be utilised. Three regions were 
selected for MOS spectroscopy to fulfill the following criteria: 
sample the inner and outer regions at a range of effective radii to 
detect any abundance gradient from the PNe; sample the major 
and minor axes; ensure that some of the brightest PNe were ob- 
served to maximize the probability of detecting the faintest line 
species; collect as many PN spectra as possible. Not all of these 
criteria are mutually consistent, but the three regions selected at 
radial offset distances of 3.4, 8.7 and 15.4' provided spectra of 
51 PNe at a rangeof m 50 07A ( m 50 o7A = -2.5 log F A (erg cirT 2 s~') 
-13.74; Jacobv fT989i ) from 23.5 to >27.0 mag. 

2.1.1. Pre-imaging 

Direct images with FORS1 and a narrow band [O III] filter 
(called OIII+50, centred at 5005A and FWHM 57A) were ob- 
tained in December 1999- January 200CQin service mode. Tab.Q] 
lists details of the imaging observations whose primary purpose 
was to provide images for the MOS slit assignment using the 
FIMS tool. The standard resolution collimator was used giving 
a field of 6.8x6.8' with a pixel size of 0.2" per pixel. The field 
names are taken from those of Hui et al. (1993a). A compan- 
ion filter (actually a redshifted [O III] filter OIII/6000, centred 
at 5109 A with FWHM 61 A) was also used to obtain an [O III] 
off-band subtraction to confirm the reality of the PNe; all the pre- 
viously catalogued PNe in the fields were confirmed as [O III] 
emission line objects. Seeing at the time of observations was be- 
tween 0.5 and 1.2". 

2.1.2. MOS Spectroscopy 

FORS1 MOS observations in service mode were conducted in 
three sessions^. Tab. [2] lists the salient details in chronologi- 
cal order. To cover the wavelength range with useful diagnos- 
tic emission lines (essentially from [O II]3726,3729A to be- 
yond [O II]7320,7330A), two grisms were employed: grism 
600B (called GRIS.600B+12) has a dispersion of 1.2A/pixel 
and, for a centred MOS slitlet, a wavelength coverage of 3450 
to 5900A; grism 300V (named GRIS_300V+10 and used with 
a GG435 blocking filter) has a lower dispersion of 2.7A/pix 
and a coverage of 4450 to 8650A. The overlap region between 
both spectra is 4500-5900A thus allowing at least the strong 
[O III]4959,5007A lines to be used to tie the two spectra to a 
common flux scale. The slit width for the campaigns in 2000 and 

1 ESO programme 64.N-02 19 

2 ESO programmes 64.B-0219, 66.B-0134 and 71.B-0134 



Table 3. FORS 1 spectrophotometric standard star observations 



Star 






Date 






(sec) 




LTT 7379 


600B 


25 


2000 05 02 


LTT 7379 


300V 


7 


2000 05 02 


LTT 7379 


600B 


25 


2001 03 24 


LTT 7379 


600B 


25 


2001 03 27 


EG 274 


600B 


30 


2003 04 08 


LTT 7379 


600B 


25 


2003 04 24 


LTT 7379 


300V 


7 


2003 04 24 


LTT 7379 


600B 


25 


2003 04 30 


LTT 7987 


600B 


100 


2003 05 05 


LTT 7987 


300V 


30 


2003 05 05 


LTT 7987 


300V 


30 


2003 06 02 


LTT 7379 


300V 


7 


2003 07 22 



2001 was 0.8" whilst the observations in 2003 had a slit width 
of 1.0". The resulting spectral resolutions are 4.8 and 6. OA for 
600B and 10.7 and 13.4A for 300V observations respectively. 
The supporting observations of the spectrophotometric standard 
stars for flux calibration, which were taken with a broad slit of 
5" width, are listed in Tab. [3] 

3. Reduction and analysis 

All the data frames (flats, science frames on PN and standard 
stars) were bias subtracted using master bias frames provided 
by the ESO reduction pipeline contemporaneous with the ob- 
serving data. Standard IRAR3 routines were used for the re- 
duction. Pixel-to-pixel flat fields were constructed from dome 
flats by extracting the 2D area of each slitlet, collapsing in the 
cross-dispersion direction, box car smoothing in the dispersion 
direction and dividing the dome flat by the smoothed version. 
Identical reductions were performed for the bluer (600B grism) 
and redder (300V grism) spectra. Wavelength calibration was 
achieved by fitting 4th order Chebyshev polynomials to the de- 
pendence of pixel position on wavelength for the arc lamp lines 
on each slitlet separately. The separate slitlet images, rebinned 
to an identical wavelength scale, were then recombined into a 
single image which was corrected for atmospheric extinction. 
Individual exposures were combined with clipping to remove 
discrepant pixels caused by detector bad pixels and cosmic ray 
events. Flux calibration was achieved by observations of one or 
more spectrophotometric standard stars, the ones applied being 
listed in Tab. [3] The standard star spectra were analysed in an 
identical way to the PN, except that the one slitlet pertaining to 
the standard star was reduced. In the few cases where a stan- 
dard star was not available in the same configuration as the PN 
on the night of observation, a standard star from another night 
had to be used. In general conditions were good or photometric, 
and no large flux calibration discrepancies were found between 
spectra of the same field observed in different runs (see Tab. [2). 
Narrow band magnitudes for the spectrophotometric standards 
were taken from Hamuy et al. (1994) and the PN spectra were 
flux calibrated using iraf.noao.onedspec routines. Fig. Q] shows 
an example of the reduced multi-slit spectra for field F56. The 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. FORS1 pre-imaging observations 



Position 


a 


S 


[O III] Exp 


Cont. 


Date 


(radius, kpc) 


h m s 


o 1 It 


(sec) 


(sec) 




F56 (9.6) 


13 25 50.6 


-43 06 09.2 


2x480 


2x220 


1999 12 27 


F42 (3.8) 


13 25 09.0 


-43 02 29.3 


2x480 


2x220 


2000 01 13 


F34(17.0) 


13 24 47.2 


-43 13 32.2 


2x480 


2x220 


2000 01 13 



Table 2. FORS 1 PN spectroscopic observations 



Position 


a 


5 


Grism 


Exp. 


Date 


Seeing 




h m s 






(sec) 




(") 




F56 


13 25 50.6 


-43 06 09.2 


600B 


4x2400 


2000 05 02 


0.8- 


1.2 


F56 


13 25 50.6 


-43 06 09.2 


300V 


2x2400 


2000 05 02 


0.8- 


1.2 


F42 


13 25 09.0 


-43 02 29.3 


600B 


2x1500 


2001 03 22 


0.8 




F34 


13 24 47.2 


-43 13 32.2 


600B 


2x1500 


2001 03 24 


0.9 




F34 


13 24 47.2 


-43 13 32.2 


600B 


2x1500 


2001 03 27 


0.8 




F42 


13 25 09.0 


-43 02 29.3 


600B 


2x1500 


2001 03 27 


0.7 




F34 


13 24 47.2 


-43 13 32.2 


300V 


2x1320 


2003 03 25 


0.7 




F34 


13 24 47.2 


-43 13 32.2 


600B 


8x1320 


2003 04 08 


0.5 




F42 


13 25 09.0 


-43 02 29.3 


300V 


2x1320 


2003 04 24 


1.0 




F42 


13 25 09.0 


-43 02 29.3 


600B 


3x1320 


2003 04 24 


1.0 




F42 


13 25 09.0 


-43 02 29.3 


600B 


3x1320 


2003 04 30 


0.6 




F42 


13 25 09.0 


-43 02 29.3 


300V 


2x1320 


2003 04 30 


0.5 




F56 


13 25 50.6 


-43 06 09.2 


300V 


2x1320 


2003 05 05 


0.5 




F56 


13 25 50.6 


-43 06 09.2 


600B 


4x1320 


2003 05 05 


0.6 




F34 


13 24 47.2 


-43 13 32.2 


300V 


2x1320 


2003 06 02 


1.1 




F56 


13 25 50.6 


-43 06 09.2 


600B 


2x1320 


2003 07 21 


1.0 




F56 


13 25 50.6 


-43 06 09.2 


300V 


2x1340 


2003 07 22 


0.4 





unequal wavelength coverage of the spectra is dependent on the 
relative position of the target centred slitlet within the field. Thus 
a target to the right edge of the field is truncated to the red but 
extends to lower wavelengths than a target in the centre of the 
FORS 1 MOS field. 

Fig. [TJ demonstrates that there can be more than one PN 
per slitlet; in some cases this was planned, by placing the slit- 
let such that two PN from the Hui et al. catalogue lay on the 
slit, but in some cases in the two inner fields, PN were spec- 
troscopically detected that were not present in the Hui et al. 
catalogue, on account of their faintness or proximity to an- 
other catalogued PN. These previously uncatalogued PN typi- 
cally had ni5oo7A > 27mag. On the basis of the strongest line - 
[O III]5007A - the number of PN detectable in the three fields 
were: 21 in F56; 21 in F42; and 9 in F34. On each slitlet, regions 
distinct from the PN spectra were located and used to subtract 
the background, with a 1 st or 2nd order fit in the cross-dispersion 
direction. In determining the spectra of the PN, no attempt was 
made to perform a separate sky background subtraction since 
there is galaxy background at all slitlet positions, except perhaps 
for the slitlets at largest galactocentric distance in the outermost 
field (F34). ID spectra of each PN were formed by summing 
the background subtracted signal along the slit; 7-9 pixels (1.4- 
1.8") were used to collect around 60 to 70% of the flux for 0.9" 
seeing (assuming a Gaussian PSF). In addition to producing ID 
spectra of flux, the same extraction was performed for non-flux 
calibrated spectra in order to determine the random noise errors 
on the extracted spectra. No attempt was made to propagate other 
sources of error, such as flat fielding, into the resultant error vec- 
tors 



As is evident from Tab. [2] all fields were observed on more 
than one occasion. In order not to lose any information, all the 
flux calibrated and extracted ID spectra were employed to form 
averaged spectra per PN. It was expected that the average should 
be formed using the inverse of the exposure time as weight; 
however the sky transmission, airmass, seeing and moon phase 
could all affect the resulting noise in a spectrum, so a simple un- 
weighted mean was finally employed. However the F56 600B 
spectra obtained in 2003 were found to have lower fluxes but 
similar signal-to-noise (S/N) as the 2000 observations. These 
later observations were rescaled in forming the average. For each 
PN, the various spectra were intercompared before averaging; if 
a line was found to be significantly discrepant between spectra, 
a careful examination of the raw 2D spectra was made to deter- 
mine if a processing step, such as incorrect CR removal in the 
PN spectrum or the local background, had affected the emission 
line flux. If a line was deemed to be badly affected, the final 
spectrum was substituted using only the clean spectrum. 

The 300V spectra overlap with the 600B spectra for the 
range 4600 - 5800A allowing the red lines (He I 5876A, 
Ha'+[N II], [S II], etc) to be placed on the same relative flux scale 
as the blue lines, primarily using the strong 5007 A line to scale 
the spectra. The errors on the final spectra were combined using 
Gaussian error propagation though these combination steps . The 
resultant spectra thus have a range of line S/N and the magni- 
tude of the propagated errors were checked in two ways: regions 
of line-free continua were chosen and the root-mean square on 
the mean was compared with the mean of the statistical errors 
for the same regions; the rms on the mean value of the fixed 
[O III]5007/4959A ratio was compared with the mean of the 
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E 
t 




Fig. 1. The reduced MOS spectra for field F56 taken with grism 
600B (data from May 2000). Eighteen of the 19 slitlets are 
shown for a wavelength range from about 3970 to 5520A and 
the numbering of the extracted PN spectra are indicated. Two 
slitlets show multiple PN spectra. The higher background stel- 
lar continuum towards the top of the image corresponds to lower 
galactocentric distance. 



errors on this observed ratio from the Gaussian line fits (see 
Section 4.1). It was found that the naively propagated errors 
over-estimated the real errors on the data values as demonstrated 
by these two tests; the factor varied slightly between data sets 
but was around 1 .7 for each spectrum. The statistical errors in 
each spectrum were amended by this amount and are those that 
are listed as the propagated errors on the measured line fluxes 
(see Tab. [IB- 

In order to examine the stellar spectra at the positions of the 
PN, an attempt was made to subtract the sky background for 
the bluer (600B) spectra. The F34 field contains the highest pro- 
portion of sky over galaxy background, so was used to remove 
sky from the exposures of the other fields when F34 and other 
fields were observed on the same night (observing runs in 2001 
and 2003 - see Tab. [2). The flux calibrated sky spectrum per 
pixel for the F34 field, excluding the PN spectra, was formed 
and subtracted from the F56 and F42 data. Some mismatches in 
terms of the sky lines were noted, particularly [N I], and some 
of the sky spectra produced poor subtraction in the sense of 
over-subtraction to short wavelengths. These problematic spec- 



tra were not used in forming a mean galactic background spec- 
trum in the two fields F42 and F56. 



4. Results 

4. 1 . Individual PN spectra 

The ID spectra of the 51 PN in NGC 5128 were analysed by 
interactively fitting Gaussians to the emission lines with a lin- 
ear interpolation to the underlying galaxy continuum over the 
line extent. Errors on the line fits were propagated from the flux 
errors. The extinction correction was determined by comparison 
with the case B values for 12000K and 5000cm 3 and the Seaton 
( 119791 ) Galactic reddening law with R=3.2 (in the absence of 
other information on the appropriate reddening law to adopt for 
NGC 5128), in all cases where at least the Ha and H/3 lines were 
detected. The observed line fluxes (normalised to Hj8=100) and 
errors are presented in Tab. [TTJ The field number and slitlet num- 
bering of the target is provided, the number from the catalogue 
of Hui et al (Q~993b) together with the Ac, A6 offsets (in arcsec) 
from the position of the nucleus (taken as the SIMBAD coordi- 
nate 13* 25 m 27.6 s -43° 01' 08.8" (J2000)), the observed loga- 
rithmic H/3 line flux and error, and the 1115007A determined from 
the observed [O III]5007A line flux and error. The extinction 
correction c, and Eb-v values are listed in Tab. [12] The errors on 
the extinction were not propagated to the dereddened line errors. 

Seven PNe with weak emission lines, whose large errors on 
the H/3 line precluded reliable determination of line ratios, were 
also analysed but are not included in Tab. [TT] these are listed 
separ ately in Tab. [4] Except for the PN F56#7 (5509 in Hui 
et al. I1993bl) . the absolute coordinates of these faint PN were 
determined from the relative positions of the PN on the slitlet 
with respect to the PN coordinates from Hui et al.((T993b'i. The 
[O III]/HjS ratios are listed in Tab. |4] where H/3 was detected; 
however given the large uncertainty on the H/3 line flux, these 
ratios are poorly determined. The PN not detected from the nar- 
row band imaging of Hui et al. (11993bt are all faint (1115007,4 > 
26.5) and are near the centre of the galaxy where the stellar con- 
tinuum is high. Slit spectroscopy allows lower equivalent width 
emission lines to be detected and can thus probe deeper than fil- 
ter imaging. 

The comparison between the 1115007A magnitudes from this 
work, converting the measured 5007 A flux to 0150074, and those 
from Hui et al. (1993b) is shown in Fig. [2] For Field F56, the 
observed 1115007A magnitudes are on average 0.33+0.12 mag. 
brighter than the Hui et al. measurements, excluding the two 
faintest PN. The absolute H/3 flux calibration was adopted from 
the May 2000 600B observations; presumably the spectrophoto- 
metric standard was taken under slightly poorer conditions , lead- 
ing to an apparent higher absolute flux for the PN observations. 
For the F42 field, the 1115007A magnitudes are 0.67+0.21 mag. 
fainter than the Hui et al. ones, which can be attributed to an 
offset of the slits from the PN positions, since all the F42 obser- 
vations display lower 1115007A than the imaging observations. For 
the F34 field, the 1115007A magnitudes are on average 0.51+0.21 
mag. fainter than the imaging magnitudes, which are consistent 
with a xQ.2" slit offset for 0.8" wide slits in 1.0" seeing (assum- 
ing a Gaussian seeing profile). 

4.2. Region-averaged PN spectra 

Single spectra consisting of the sum of all the observed brighter 
PN (i.e. those contained in Tab. [TT1) in each of the three fields 
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Table 4. NGC 5128 PN without reliably measured H/3 fluxes 



MOS Slitlet Hui et al. 


a 


6 


Offsets 


m5007A 


[O III]/H/? 


No. 


(h m s) 


(° ' ") 


Aa A6 


mag.* 


approx. 








(arcsec) 






F56#7 5509 


13 25 56.41 


-43 06 23.2 


(+293.6,-314.4) 


26.2 1 " 


29 


F56#13a 


13 25 46.50 


-43 04 57.3 


(+207.3,-228.5) 


27.6 


4 


F56#13c 


13 25 45.78 


-43 04 53.3 


(+199.4,-224.5) 


27.5 


13 


F42#12a 


13 25 15.21 


-43 03 12.5 


(-135.9,-123.7) 


26.5 


23 


F42#14a 


13 25 18.20 


-43 02 23.1 


(-103.1,-74.3) 


26.8 


22 


F42#14c 


13 25 19.24 


-43 07 19.9 


(-91.7,-71.1) 


27.1 


9 


F42#16b 


13 25 23.13 


-43 02 49.6 


(-49.0,-100.8) 


27.2 


9 



Notes. * Spectroscopic magnitude, corrected for slit/pointing losses. 
'•' Hui et al. 1993b]list m 50 07A of 25.76. 



t 1 1 1 1 1 1 1 1 1 1 1 1 r 




Fig. 2. Comparison of the imaging m^ooiA magnitude from Hui 
et al. ( 1993b) with the equivalent magnitude computed from the 
observed [O III]5007A spectroscopic flux converted to magni- 
tude. The PN in the three fields are shown with different sym- 
bols, as indicated. 



were constructed. The primary motivation was to enable the de- 
tection of fainter diagnostic lines (such as He I, [S II], [Ar IV], 
etc.) which were not detected, or were marginally detected, on 
the individual spectra. This was also an explorative study to mea- 
sure how different the derived abundances would be from the 
summed spectra in comparison with the ensemble of the indi- 
vidual PN abundances. This aspect is of interest for studies of 
more distant PN populations where all line detections from indi- 
vidual PN spectra have low S/N and summed spectra are manda- 
tory for measurement of abundances. The radial velocity of each 



PN was measured from the strong lines in the spectrum (princi- 
pally [O III] 4959 and 5007 A and Ha) and employed to shift all 
the spectra in a single region to zero radial velocity. The spectra 
were then summed and the resulting emission line spectra fit- 
ted in the same way as the individual PN spectra (Section 4.1). 
Figure 2 of Walsh et al. (2005) shows the resulting three spec- 
tra. Tab.|5]lists the resulting observed fluxes. The extinction was 
calculated from the H line ratios with the same assumptions as 
for the individual PN spectra (Sect. 4.1) and the dereddened line 
fluxes are listed in Tab. [6] 

The ionic abundances of the elements corresponding to the 
detected species can be determined from the dereddened line 
ratios with respect to H/3 (Tab. [SJ. The electron temperature 
of the ++ emitting region can be directly measured from the 
5007/4363A, but also required for abundances is an electron 
density measure. The only directly determined N e value is from 
the [S II]67 16/6731 A ratio, which has two disadvantages in 
that it samples the minority low ionization nebular volume and 
has a relatively low critical density for collisional de-excitation. 
Coppeti & Writzel (1989) show that [S II] N e measured in PN 
is similar to N e measured for higher ionization species, such as 
[CI III] and [Ar IV]. Adopting a value of N e of 5000 ctrT 3 for 
the NGC 5128 PN appears to be fair and will not lead to bias on 
the abundance estimates, unless the density is very high (^20000 
ctrT 3 ). In the high density case, the [O III]5007/4363A line ratio 
is decreased by collisional de-excitation, leading to too high a 
value of T e . 

Tab.|7]lists the derived ionic abundances for the summed re- 
gion spectra. The sources for the atomic data for the collisionally 
excited lines were taken from the compilation in Liu et al. (2000) 
and the routines from Storey & Hummer (1995) were used for 
the recombination line emissivities of H and He. The correc- 
tions for unseen stages of ionization (ionization correction fac- 
tors, ICFs) were taken from Kingsburgh & Barlow (1994). The 
errors in the total abundances do not take into account the errors 
in the electron temperature; if these are propagated the errors 
rise to +0.10 for the O abundance. The He/H abundance is only 
listed if the lines at 4471 and or 5876A were detected; if He II 
4686A was detected in addition, it was added to derive the total 
He/H. If only He II was detected no He/H is listed. 

4.3. Stellar absorption lines 

The slitlets sample a spectroscopic background consisting of sky 
with a substantial contribution from the galaxy continuum, ex- 
cept for the outermost field, F34. The regions of the slitlets not 
occupied by planetary nebula spectra can provide the spectrum 
of the stellar continuum of NGC 5 128. The inner fields, F56 and 
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Table 5. Summed PN spectra of three fields in NGC 5128 - Observed line fluxes 





F56 




F4z 




F34 




Species 


A (A) 


Fobs 


± 


Fobs 


± 


Fobs 


+ 


[Oil] 


3727 


28.6 


1.7 


23.9 


2.2 


54.5 


2.9 


[Ne III] 


3868 


84.0 


1.7 


62.3 


1.5 


74.0 


1.9 


HI 


3889 


17.3 


4.4 


13.3 


1.6 


14.3 


1.2 


[Ne III] + He 


3970 


33.5 


1.0 


26.9 


1.5 


34.4 


1.9 


H<5 


4101 


21.1 


1.2 


10.7 


1.2 


21.3 


1.1 


Hy 


4340 


46.7 


1.9 


41.5 


2.0 


40.8 


1.5 


[O III] 


4363 


8.7 


0.7 


6.1 


1.3 


8.5 


0.6 


He I 


4471 


8.0 


1.0 


6.5 


1.0 


5.0 


3.3 


He II 


4686 


10.7 


1.2 


9.3 


1.9 


13.1 


1.1 


H/? 


4861 


100.0 


0.0 


100.0 


0.0 


100.0 


0.0 


[O III] 


4959 


391.3 


5.7 


396.4 


7.1 


417.1 


7.1 


[O III] 


5007 


1159.6 


16.0 


1213.6 


20.2 


1247.8 


20.6 


He I 


5016 


3.7 


0.4 










[N I] 


5199 


3.2 


0.2 










He I 


5876 


14.5 


1.9 


29.9 


5.6 






[Nil] 


6548 


50.3 


1.2 


62.7 


4.4 


62.7 


4.6 


Ha 


6562 


390.6 


3.8 


422.7 


9.8 


412.8 


14.7 


[Nil] 


6583 


152.2 


1 o 

1.8 


1 AC\ 1 

149.2 


A 1 

4.3 


1 If H 

115.7 


5.5 


He I 


6678 


4.6 


0.8 


7.1 


4.6 






[S II] 


6716 


4.5 


0.4 


7.3 


2.8 






[S II] 


6730 


10.9 


1.1 


27.9 


3.0 






[Ar III] 


7133 


28.5 


1.1 


22.9 


3.1 






[Oil] 


7325 


27.2 


2.4 


19.6 


4.6 






log F(H0) 




-15.46 


0.01 


-15.42 


0.01 


-15.78 


0.01 


ni5007A 




22.24 


0.03 


22.10 


0.04 


22.98 


0.04 



Table 6. NGC 5128 PN Region spectra - Dereddened line fluxes 





F56 




F42 




F34 




Species 


A (A) 


Fobs 


± 


Fobs 


± 


Fobs 


+ 


[Oil] 


3727 


37.0 


2.3 


33.0 


3.0 


73.8 


3.9 


[Ne III] 


3868 


105.9 


2.1 


83.3 


2.1 


97.2 


2.5 


HI 


3889 


72.0 


5.5 


17.7 


2.1 


18.7 


1.5 


[Ne III] + He 


3970 


41.4 


1.2 


35.0 


1.9 


44.2 


2.4 


H<5 


4101 


25.4 


1.4 


13.5 


1.4 


26.4 


1.4 


Hy 


4340 


53.1 


1.8 


48.7 


2.4 


47.4 


1.7 


[O III] 


4363 


9.8 


0.8 


6.6 


1.3 


9.9 


0.6 


He I 


4471 


8.9 


1.0 


7.4 


1.2 


5.6 


3.6 


He II 


4686 


11.2 


1.3 


9.8 


2.0 


13.8 


1.1 


H/3 


4861 


100.0 


0.0 


100.0 


0.0 


100.0 


0.0 


[O III] 


4959 


381.8 


5.6 


384.6 


6.9 


405.5 


6.9 


[O III] 


5007 


1118.3 


15.4 


1160.0 


20.0 


1196.3 


19.7 


He I 


5016 


3.6 


0.4 










[N I] 


5199 


3.0 


0.2 










He I 


5876 


11.6 


1.5 


22.8 


4.3 






[Nil] 


6548 


36.5 


0.9 


42.1 


2.9 


43.0 


3.1 


Ha 


6562 


282.9 


2.8 


282.8 


6.6 


282.8 


10.1 


[Nil] 


6583 


109.9 


1.3 


99.4 


2.9 


79.0 


3.7 


He I 


6678 


3.3 


0.5 


4.7 


3.0 






[S II] 


6716 


3.2 


0.3 


4.8 


1.8 






[S II] 


6730 


7.7 


0.8 


18.1 


1.9 






[Ar III] 


7133 


19.1 


0.8 


13.9 


1.9 






[Oil] 


7325 


17.9 


1.6 


11.6 


2.7 






c 




0.44 


0.03 


0.55 


0.07 


0.51 


0.12 


log F(Hj8) 




-15.03 


0.01 


-14.88 


0.01 


-15.27 


0.01 


ni5007A 




21.18 


0.03 


20.74 


0.04 


21.74 


0.04 



F42, contain no slitlets sampling the sky free of stellar contin- 
uum, and so no simultaneous sky subtraction can be performed. 
For the outer field F34 however, the stellar continuum is very 
weak and was not detected, as shown by the absence of a gra- 
dient in the background with radial offset from the galaxy cen- 
tre. The 600B F34 spectrum taken on 2003-04-08 was thus used 



to produce a candidate mean sky spectrum for subtraction from 
the spectra for the two inner fields (F56 and F42). This sky 
was interactively subtracted from the F42 and F56 spectra for 
each slitlet. Good results were found only for the F42 spectra 
taken in April 2001 and for the F56 field spectra taken in May 
2003. These results are attributable to very differing sky contri- 
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Table 7. NGC 5128 PN Region spectra - Diagnostics and Abundances 



Species 


EC/; 

rjo 




17 A O 

r4Z 




~C2 A 

ro4 






Value 


+ 


Value 


+ 


Value 


+ 


[O I11J (5(X)7+4959)/4363A 


153.1 


12.6 


234.0 


46.2 


161.8 


10.0 


I ( \ TTTI T / T/- \ 

|U 111J i e (K) 




300 




560 


lU/oU 


Tin 
zzU 


(J 


1. / Oil- J 


1 IE t 

l.lE-o 


I.oyh-D 


2.4E-6 


3.jze-D 


l.ob-0 


LNC 


7 Q^F S 


1 f\ 
1 .uc u 


1 fHF 4 




7 9SF S 


1 QF f\ 


FT^+ 


n 1 1'x 

V. 1 1 J 








1 HQ 


n n7n 

U.U l\) 


He (4471 A) 


0.0094 


0.0011 


0.0081 


0.0016 


0.0115 


r\ f\r\f\c\ 

0.0009 


o ++ 


3.091E-4 


4.3E-6 


4.997E-4 


8.6E-6 


3.307E-4 


5.4E-6 


He + (5 876 A) 


0.084 


0.011 


0.160 


0.031 






N + 


1.83E-5 


2.2E-7 


2.28E-5 


6.6E-7 


1.32E-5 


6.2E-7 


c + 

a 


4.40E-7 


4.6E-8 


1.39E-6 


1.5E-7 






Ar ++ 


1.50E-6 


6.3E-8 


1.46E-6 


2.0E-7 






12 + Log(0/H) 


8.53 


0.01 


8.74 


0.01 


8.59 


0.02 


He/H 


0.182 


0.019 


0.150 


0.023 


0.122 


0.0070 


Ne/O 


0.257 


0.006 


0.206 


0.006 


0.220 


0.007 


N/O 


1.04 


0.07 


0.85 


0.08 


0.38 


0.03 


12 + log(S/H) 


6.55 


0.05 


7.06 


0.04 






12 + log (Ar/H) 


6.49 


0.04 


6.44 


0.06 







butions to the galaxy spectra at the different observation epochs. 
For the F56 field, the sky-subtracted data show very low con- 
tinuum in the five slitlets (1-5) furthest from the galaxy centre; 
these slitlets were not considered for analysis of the galaxian 
stellar continuum. Fig. [3] shows the resulting stellar continuum 
from fields F56 (blue line) and F42 (red line), for all the slitlets 
summed. Strong Call H & K, Hy, Hfi and Mg II lines typical of 
an intermediate-old stellar population are clearly visible and are 
indicated on the Figure. 

The absorption lines in the individual galaxy continuum 
spectra for each slitlet were analysed by computing Lick in- 
dices (Worthey et al. |1994| l. The spectra were smoothed with a 
Gaussian to simulate the ~8A resolution of the Lick/IDS stel- 
lar spectra and shifted to zero velocity before the equivalent 
widths were computed. Errors on the indices were computed by 
propagating the flux errors in the equivalent width determina- 
tion. A single radial velocity of 580 kms -1 was used to shift 
the spectra to rest wavelength. Errors of +50 kms~' in this ra- 
dial velocity cause an offset in the bounds of the Lick indices 
and can introduce errors up to about 0.2A in EW. Three of the 
Lick indices (H/3, Mgb and the combined Fe index <Fe>, de- 
fined as 0.5*[EW(Fe 527 oA) + EW(Fe 5335A )], are plotted in Fig.g] 
as a function of the projected position of the slitlet centres from 
the centre of the galaxy. The <Fe> index shows no radial gra- 
dient nor does the H/3 index, although the latter shows higher 
scatter to larger projected distance from the nucleus. The Mg 
b index shows a weak negative radial gradient, typical of early 
type galaxies (c.f. Davies et al. l 1995b . 



5. Chemical abundance determination 

5. 1 . Extragalactic PN abundances 

The electron temperature (T c ) is required In order to determine 
reliable chemical abundances for extra-galactic PNe from the 
collisionally excited lines on account of the exponential depen- 
dence of line emissivity with temperature. The electron density 
(N e ) is also required since these lines have collision cross sec- 
tions dependent on density. In order to avoid large corrections 
for unseen ionization species a range of lines of different ioniza- 
tion from neutral and singly ionized up to high ionization is also 





Ca n 


hy 


Hp Mg b " 




H&K 








H 


i .... i 


i 



4000 4500 5000 5500 

Wavelength (A) 



Fig. 3. The observed spectra of the galaxy stellar continuum, 
summed along the whole slit but excluding the positions occu- 
pied by the PN, are shown for fields F42 (upper spectrum, in red) 
and F56 (lower spectrum,in blue). The strongest absorption lines 
are identified. 

desirable; however the well-established empirical technique us- 
ing ICFs can compensate for the lack of some ionization species 
from the optical wavelength range. As is typical for faint and dis- 
tant PNe such as observed here, these criteria are not fully met. 
In the absence of measurements of weak T e and N e diagnostic 
line ratios (such as [O III] 5007/4363A and [S II]67 16/6731 A 
or [Ar IV]471 1/4740A), the practice adopted by Stasiriska et al. 
( 1998) was to use the upper limits to the strength of these faint 
lines to constrain T e and N e within reasonable limits, taking typ- 
ical values for PN from studies in the Milky Way and Magellanic 
Cloud s (e.g. T e -12000K, N e -5000 cm 3 ). Jacoby & Ciardullo 
( 119991) introduced a different approach to abundance determina- 
tion for extra-galactic PNe (in this case for M31): that of pho- 
toionization modelling using the very well-established Cloudy 
(Ferland et al. 1998) code. The gains of this latter method are 
that both strong lines with low errors, weak lines with large er- 
rors and upper limits can all be used in arriving at a satisfactory 
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Fig. 4. The Lick indices, EW(H/3), EW(Mgb) and <Fe> (defined 
as 0.5*(EW(5270A) + EW(5335A)), determined from the stel- 
lar continuum in each slitlet, are plotted as a function of the 
projected distance from the nucleus of NGC 5128, in arcsec 
(lower axis) and kpc (upper axis). Measurements for field F42 
are shown in red and for field F56 in blue. 



model, and the parameters of the central stars (luminosity, £, and 
effective temperature, T e ff) are derived as part of the photoion- 
ization modelling. 

Extra-galactic PNe are spatially unresolved but have known 
distances and absolute fluxes available though m 50()7 ^ photome- 
try. An input data set for photoionization modelling can be as- 
sembled, given some simplifying assumptions in the absence of 
more detailed information; in particular, the nebular geometry 
and information on the stellar atmosphere are not available. A 
black body has been shown to be an acceptable assumption for 
PN central stars (Howard et al. 1 19971 . which leaves the nebular 
geometry. Taking a spherical shell of constant density as a base- 
line, Jacoby & Ciardullo (1 19991 1. hereafter JC99, compared neb- 
ular abundances using directly measured T e and N f values and 
ICFs to check that photoionization modelling was indeed a valid 
approach. The approach was shown to work well. Additional 
complexity, such as two zone density structure and the presence 
of dust within the nebula, was found to be required to satisfacto- 
rily match the observed spectra in some cases. 

Magrini et al. (2004) followed the same procedure as JC99 
in determining abundances of PNe in M33. They also tested the 
method by applying it to the integrated spectra of well-observed 
Milky Way PNe and compared the model abundances to those 
from the ICF method, finding agreement within 0. 15 dex for O/H 
and 0.3 dex for N/H. In addition Magrini et al. (j2004) compared 
abundances derived from Cloudy models with those from ICFs 
for three PNe of their M33 sample. They also considered com- 
parisons employing model atmospheres and an r~ 2 density in the 
nebular shell, and found no major discrepancies. The technique, 
involving the development of a large number of models in work- 
ing to a satisfactory match of observed and model spectrum, is 
labour-intensive and therefore not suitable for large samples. It 
also requires simplifying assumptions. The technique, however, 
is well suited to the present sample of PN spectra in NGC 5128 



where the range of abundances is of particular interest and the 
data set is limited; it also serves as a probe of the star formation 
history. 

Only 10 of the PN spectra presented in Tables QT| and [12] 
have detectable [O III]4363A, enabling direct measurement of 
the electron temperature, and only 7 have an electron density 
sensitive ratio ([S II]6716/6731A) measured. The typical S/N 
on the weak 4363A line (< 3) implies that T e errors of 1250K 
at 11000K affect abundance determination, leading to errors of 
on log (0 ++ /H + ). For the PNe without detectable diagnos- 
tic line ratios, two approaches are possible: employ the region- 
averaged spectra to assign T e and N c values to the single PN, 
although these values may be substantially in error for particu- 
lar nebulae; or, employ photoionization modelling of the strong 
lines to find the best-matching model that does not violate the 
upper limit constraints for the weak lines. 



5.2. Modelling procedure 

The simplest possible photoionization model was employed: a 
black body central star characterized by its temperature and lu- 
minosity, the latter given by the H/3 luminosity derived from 
the 1115007A magnitude (Hui et al. 1993b), the reddening and the 
dereddened 5007A/HjS; a spherical shell with the inner radius 
set to a small value (0.005 pc) to ensure high ionization emis- 
sion close to the central star, but not so close that instabilities 
arise in the model; the outer radius set to a large value (0.5 pc) 
since the nebulae are expected to be optically thick, as given the 
high luminosity of the PNe observed in NGC 5128 (i.e. within 2 
mag. of the peak of the PNLF); and an initial set of abundances. 
For the abundances of He, N, O, Ne, S and Ar, the initial values 
were taken from the ICF analysis for the integrated spectra listed 
in Tab. [8] and the C/O ratio taken as 0.5. However these values 
were immediately allowed to vary to fit the individual spectra 
(with C varying in lock-step with O) and should not be seen as 
prejudicing the individual determinations from the higher S/N 
integrated spectra. The initial assumption of an optically thick 
nebula could also have been relaxed in the modelling process 
but was not found to be demanded by the fits to the PN spectra. 
The philosophy was to aim to match the dereddened lines fluxes 
within the listed errors for all the spectra listed in Tab. [T2~l(except 
F56#9, F56#17, F42#15b and F42#17 without detected [O II] or 
He lines), thus 40 PN spectra. These PNe represent the upper 2. 1 
mag. of the PNLF. 

The modelling procedure followed closely that outlined by 
JC99 and also Marigo et al. 2004 An initial default density of 
5000cm -3 was used. The initial estimate of T e ff was taken from 
He II 4686A/Hj6; in the absence of 4686A, its upper limit and 
the 5007A/H/? ratio was used. The presence of 4686A enabled 
T e ff to be refined within a few runs of Cloudy; if 4686A was not 
present, both T e ff and O/H had to be adjusted. Then He was ad- 
justed to match the He I lines, with minor adjustments to T e for 
the concomitant changes in He II flux. O/H and the density were 
adjusted to match [O II]3727A and, if present, the [O IIJ7325A 
quadruplet. For large changes in O abundance, Ne and C were 
initially changed in proportion. Once the basic physical param- 
eters (T e ff, density, He, O) were roughly determined, Ne, N, Ar 
and S were altered to match the observed lines. In the final phase 
all lines and T e ff and density were subject to small changes to 
improve the fit, where possible. A guide to the goodness of fit 
of the observed and predicted line strengths was formed, taking 
the sum of the line flux differences normalised by the S/N: thus 
larger discrepancies for weak lines could be weighted lower. A 
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discussion of the goodness of the matches and the estimation of 
errors in presented in section 5.4. In general around 10 iterations 
per PN was enough to reach a satisfactory match to the spectrum 
but in more difficult cases (F56#12b, F42#10 and F34#7 for ex- 
ample) up to 30 iterations were required. In total, ~600 Cloudy 
models were run to complete this analysis. 

5.3. Results and scrutiny 

Tab.[8]lists the derived abundances for He, N, O, Ne, Ar, and S 
with respect to H where lines of these elements were detected, 
together with the stellar luminosity, black body temperature and 
the constant shell density. The quality of the derived parameters 
are broadly classified into three categories based on: He II (and 
or [O III]4363A) and [O II] both detected (15 PN: grade a); [O II] 
but not He II detected (18 PN: grade b); neither He II nor [O II], 
but at least He I, [O III]5007A, [Ne III]3868A and [N II]6583A, 
detected (6 PN: grade c). Ten of the quality "a" spectra had 
4363A detected and the models tailored to fit the strong lines 
(i.e. a model was not tweaked to exactly fit the 4363A line) were 
found to match the 4363 A line within the lcr errors in all cases 
except one (F42#16a; 4207 in Hui et al. U993bb . The results pre- 
sented are not claimed to be unique, as the adopted geometry 
of the nebulae is the simplest possible and the use of a black 
body atmosphere does not produce the most satisfactory fits to 
the level of ionization in detailed photoionization modelling of 
Galactic PNe (e.g. Wright et al. 1201 11 1. However these simple 
models are capable of matching the observed variety of spectra 
with a very plausible range of abundances, central star tempera- 
tures and nebular densities. 

Carbon is the most important coolant after O, but no lines 
of C were detected so there are no direct constraints on the C/H 
abundance. In the absence of other evidence, the abundance of C 
was assumed to be similar to O (C/O=0.5) and to vary in lock- 
step with O as the abundance of O was altered to match the spec- 
trum. Once an adequate model had been found, tests were per- 
formed on four of the spectra, altering the C abundance by a fac- 
tor 2.5 and refitting for the other species. A maximum difference 
in 12+log(0/H) of -0.10 was found for 12+log(C/H) higher by 
+0.40; for 12+log(C/H) lower by 0.40, a maximum difference 
of +0.02 in 12+log(0/H) was found. Smaller changes in N and 
Ne abundances were required to match the spectra with these al- 
tered C abundances. Thus the abundances are not very sensitive 
to modest changes in the assumed C/O ratio. 

For the best observed PN spectrum (viz. the target with the 
brightest m 50O7A ) PN, (F56#2; 5601 in Hui et al. U993bl) . a set of 
Cloudy models were run with NLTE model atmospheres from 
the Tubingen compilation (Rauch 2003), and also with the pres- 
ence of dust inside the ionized region. H, He and PG1 159 model 
atmospheres were employed. While the H model atmosphere 
showed a much lower T e ff, an He or PG1 159 atmosphere at 170- 
180 kK produced satisfactory fits, with the resulting O/H lower 
by up to 0.16 dex. Including dust with an atmosphere model re- 
quired that O/H be reduced by up to 0.05 for a gas-to-dust mass 
ratio of 1000; as the primary coolant of the nebula, the lower 
O/H is compensated by the enhanced cooling from dust. 

Stasihska (2002) questioned the uniqueness of abundance 
determinations even in the case that a direct measure of T e from 
e.g. [O BTJ4363A is available. Depending on the temperature of 
the central star, a strong gradient in the electron temperature can 
exist within the nebula and a single value of T e may not be a 
good representation for the bulk of the emission. The T e gradi- 
ent is far less steep for PNe with hot central stars (~100000K) 



than for cooler central stars and the use of a single T e should 
yield reliable abundances. Stasinska (2002) presented photoion- 
ization models for a Milky Way Bulge PN with a central star 
temperature of 39 000K that can be equally well fit by a lower 
abundance ([O/H]=-0.31) and by a higher abundance model 
([O/H]=+0.39fl 95% of the PNe observed in NGC 5128 appear 
to have T e ff > 60kK, so the likelihood of a double-valued abun- 
dance is very low. However since the lowest T e ff values derived 
in the Cloudy modelling (Tab. 7) are ~50kK, then uniqueness 
may be a more significant concern for these few objects. 

A concerted attempt was made to fit the spectra of the PNe 
with the lowest derived T e ff by low and high metallicity (Z) 
models. Here there is the liberty to change the stellar temper- 
ature since the He II 4686/H/? ratio can be taken as an upper 
limit based on the errors on local lines. Two PN were found 
that could not be distinguished in terms of lower and higher Z 
Cloudy models. In the case of F56#12b (a class b spectrum with 
the [O II] line detected but not He II), the higher Z solution with 
12+log(0/H) of 8.55 for a BB T eff of 49 kK could not be dis- 
tinguished from a model with T e /f of 73 kK and 12+log(0/H) 
of 8.06 if the density was decreased from 4.0 X 10 4 to 1.8 x 10 4 
cirT 3 . For F42#3 (class c spectrum without [O II] and He II lines 
detected), a higher Z solution with 12+log(0/H) of 8.44 for T eff 
of 50 kK was indistinguishable from that with 12+log(0/H) of 
8.01 but T e ff of 73kK, for the same density in the shell. For 
the class c PN, an attempt was made to match each spectrum 
by a forced lower Z model but no convincing matches could be 
found; nevertheless, one may exist given the rather weak con- 
straints implied by these lower quality spectra. Conservatively, 
the higher abundance solutions were adopted (in Tab. [8) on the 
grounds that they do not differ strongly from the values for other 
PNe at similar radii. 

The two PNe with alternate low O/H values are not improba- 
ble and are comparable to the abundances of PNe in dwarf galax- 
ies (e.g. Dopita & Meatheringham, |T99 1 l and Leisy & Dennefeld, 
2006); F42#4 also shows a similarly low value of O/H but has 
a well-detected spectrum. In the context of a giant elliptical 
galaxy, such low abundance values are surprizing because they 
imply low metallicity progenitor stars,. Radial velocities of these 
PNe (available for F56#12b and F42#4 from Peng et al. l2004b do 
not show 'peculiar' velocities, being within lcr of the mean ve- 
locity of NGC 5128. There is no indication of these PNe being 
obvious 'halo' objects or of resulting from a recent low-mass 
galaxy interloper not yet in dynamical equilibrium with the par- 
ent galaxy. 

The abundances of the individual PN from Tab. [8] were com- 
bined and compared to the abundances derived from the inte- 
grated region spectra with empirically determined temperatures 
and ICFs presented in Tab. [7] The means were determined by 
weighting the individual abundances by rri5007A. For the com- 
bined region F56, the weighted mean abundances from Tab. [8] 
for O agreed to within 0.01 dex and for F34 to within 0.06 dex 
(weighted means of the log and the fractional O/H abundances 
were calculated with similar results). For field F42, the discrep- 
ancy on the weighted mean of the log O abundance was larger, 

4 Attempts to reproduce the photoionization fits of Stasinska ( 2002 ) 
using Cloudy 08.01 were not entirely successful. The lower metallic- 
ity model could be satisfactorily reproduced with a slightly higher T e ff 
of 43 kK and lower O of [O/H] =-1.02; however the higher metallicity 
model required [O/H]=+0.86 and at such high abundance the Cloudy 
models were very sensitive to extremely small changes in abundances 
of O, Ne, N and of the nebular size and density. From this numerical 
experiment, it is tentatively concluded that such extreme conditions are 
rare. 
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Table 8. Cloudy photoionization modelling of 40 PN in NGC 5128 



PN# 




Elemental Abundances 12 + log(A/H) 






LogL 


Ne 


Quality 1 




(He/H) 


(N/H) 


(O/H) 


(Ne/H) 


(S/H) 1 


(Ar/H) 2 


(kK) 


(£g) 


(cm- 3 ) 




F56#l 


11.06 


8.20 


8.40 


7.70 


6.55 


6.00 


60 


3.68 


5000 


b 


F56#2 


11.15 


7.85 


8.43 


7.72 




6.00 


155 


4.20 


10000 


a 


F56#3 


11.06 


8.30 


8.48 


7.78 


6.55 


6.00 


64 


3.66 


2000 


b 


F56#4 


11.06 


8.30 


8.38 


7.70 


6.55 


6.50 


78 


3.75 


5000 


b 


F56#5 


11.06 


7.40 


8.55 


7.85 




6.35 


150 


3.74 


10000 


c 


F56#6 


11.20 


8.30 


8.64 


7.95 




6.60 


75 


4.15 


6000 


b 


F56#8 


11.20 


8.62 


8.60 


7.80 


6.80 


6.50 


63 


3.77 


9000 


b 


F56#10 


11.10 


8.20 


8.38 


7.60 


6.55 


6.60 


88 


4.03 


15000 


b 


r JXJTT 1 1 


1 n on 


8.42 


8.54 


7.70 


7.00 




Qfl 


J.U / 




a 


F56#12a 


11.10 


8.53 


8.80 


8.10 


6.55 


6.40 


200 


3.86 


10000 


a 


F56#12b 


11.00 


8.80 


8.55 


7.85 




6.50 


45 


3.93 


35000 


b 


F56#13b 


11.00 


7.60 


8.52 


7.70 




6.20 


105 


3.76 


20000 


a 


F56#14 


11.06 


8.40 


8.46 


7.70 


6.55 


6.30 


95 


3.86 


15000 


b 


F56#15 


11.06 


7.90 


8.47 


7.70 


6.55 


6.20 


70 


3.86 


10000 


c 


F56#16 


11.06 


8.20 


8.47 


7.60 


6.55 


6.20 


70 


3.91 


10000 


b 


F56#18 


11.06 


8.63 


8.51 


7.85 


6.90 




90 


3.87 


15000 


b 


F42#l 


11.12 


7.85 


8.45 


7.65 


6.90 


6.20 


80 


3.83 


10000 


c 


F42#2 


11.06 


8.19 


8.68 


8.03 






155 


4.05 


10000 


a 


F42#3 


11.06 


8.26 


8.44 


7.50 




6.40 


50 


3.99 


10000 


c 


F42#4 


11.01 


8.12 


8.03 


7.10 


6.40 


6.20 


86 


4.00 


12000 


a 


F42#6 


11.15 


7.30 


8.44 


7.60 


7.10 




63 


3.63 


13000 


b 


F42#7 


11.06 


8.00 


8.51 


7.70 






127 


3.64 


10000 


b 


F42#8 


11.06 


8.20 


8.46 


7.70 






93 


3.86 


15000 


b 






8.20 


8.45 


7.70 




f\ ^n 


74. 


J. / o 


unnnn 

JwUU 


K 

U 


F42#10 


11.28 


8.18 


8.10 


7.40 


6.85 


6.10 


120 


4.29 


2500 


a 


F42#l 1 


11.06 


8.25 


8.28 


7.35 


6.85 


6.20 


100 


3.72 


10000 


c 


F42#12b 


11.06 


8.22 


8.69 


8.07 






108 


4.09 


10000 


b 


F42#13 


11.06 


8.28 


8.90 


8.25 






190 


3.84 


16000 


a 


F42#14b 


11.06 


8.20 


9.00 


8.40 


6.70 




220 


4.06 


20000 


a 


F42#16a 


11.06 


8.28 


8.72 


7.92 






100 


3.78 


14000 


a 


F42-18 


11.06 


7.68 


8.49 


7.45 






120 


4.35 


10000 


c 


F34#l 


11.06 


8.25 


8.73 


8.00 




5.65 


71 


3.89 


8500 


a 


F34#2 


11.23 


7.62 


8.09 


7.38 






90 


4.11 


1500 


a 


F34#4 


11.06 


8.10 


8.46 


7.70 




6.30 


67 


3.90 


8000 


b 


F34#7 


11.06 


8.30 


8.91 


8.30 




6.50 


82 


3.99 


13000 


a 


F34#l 1 


11.16 


7.40 


8.53 


7.78 






180 


3.59 


30000 


a 


F34#12 


11.16 


7.35 


8.48 


7.30 






167 


3.13 


30000 


a 


F34#14 


11.06 


7.90 


8.62 


7.82 






70 


3.61 


20000 


b 


F34#15 


11.06 


7.40 


8.43 


7.76 






123 


3.90 


10000 


b 


F34#16 


11.06 


8.28 


8.63 


7.95 


6.55 




180 


3.70 


5000 


a 



Notes. 1 12 + log(S/H) listed if default value of 6.55 not applied. 

2 12 + log(Ar/H) listed if default value of 6.00 not applied. 

3 Quality, a: He II and [O II] both detected; b: [O II] but not He II detected, or He II detected but not [O II]; c: neither He II nor [O II] detected. 



with a value up to 0.22 dex lower than the value in Tab. [7] The 
difference for F42 is quite striking in this comparison; the elec- 
tron temperature was measured to be ~1 000K lower than for the 
other two regions. It is suggested that the value of T f was under- 
estimated in this merged spectrum due to a poorly determined 
[O III]4363A flux. 

A further check on the integrated region spectra was per- 
formed by running photoionization models for the region spec- 
tra presented in Tab. [6] The abundances were compared to the 
empirical abundances derived using ICFs (Tab. |7). Satisfactory 
Cloudy models could not be fitted without departing from the 
simplifying assumptions used for the individual PN models. In 
particular, when the stellar temperature was chosen from the 
He HfH/3 ratio, this could not be made compatible with the elec- 
tron temperature derived from the [O III]5007/4363A ratio. If 
He II4686A matched the observations within the errors, then 



[O III]4363A was predicted much too strong (higher T c ); when 
[O III]4363A matched the model, He II was modelled too low 
(lower T e ff). The addition of modest amounts of dust inside 
the ionized region of the photoionization model to provide ex- 
tra heating through photoelectric emission from the grains did 
not resolve this discrepancy. However allowing for these mod- 
elling issues, the range of O/H was close (within 0.1 dex) to the 
empirical values listed in Tab.|7]for the sum of PNe in fields F56 
and F34, but only within 0.2 dex for F42. These demonstrations 
show that summing PN spectra provides a useful indicator of 
the mean O abundance in cases when the line detections have 
low S/N, such as for PNe in more distant galaxies. This conclu- 
sion is subject to the condition that the PN abundances do not 
differ strongly, as in this case. Mendez et al. d2005 ) have pre- 
sented lower limits on O and Ne abundances based on just such 
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a summed spectrum for 14 PNe in NGC 4697 (an elliptical at 1 1 
Mpc). 

5.4. Towards a comparative study of models 

A number of Cloudy photoionization models were run in order 
to arrive at an adopted model for each PN matched in Tab. [8] The 
models presented cannot be claimed to be unique but represent 
a feasible match to the spectra in terms of the abundances and 
physical conditions within the restricted set of parameters (black 
body atmosphere for the stars, linear density law, etc). Magrini 
et al. (120041) and JC99 discussed the accuracy of photoionization 
models to their PN spectra. In general the "a" quality spectra of 
the NGC 5128 PN approach those of the fainter sample of the 
much closer PNe observed in M33 and M31 by these authors. 
Simply adopting the range of models within the statistical errors 
on the line ratios underestimates the real uncertainties on the 
abundance determinations. For a selection of models, indepen- 
dent Cloudy runs were performed by two of the authors and the 
results compared (c.f. Jacoby et al. 1 19971 1 and with the ICF deter- 
minations in the cases where [O III]4363A was detected (using 
an adopted N e value of 5000 cm -3 ). While a full error analysis 
is outside the scope of this work, a limited investigation was at- 
tempted. A range of parameter values around the adopted ones 
were explored to determine the sensitivity of the models, given 
the errors on the spectra. 

It would be very time consuming to perform an investigation 
of the likely parameter range for all PN in Tab. [8] so two PN were 
selected. F56#2 was chosen as a high S/N case with both [O III] 
4363A and He II 4686A well detected (a class "a" spectrum in 
Table [8]l and a high stellar temperature. F34#l was chosen as a 
PN with lower S/N and only a marginal detection of He II 4686 A 
(class "b" in Table [8]) and a lower stellar temperature and repre- 
senting the typical average quality of spectra. Cloudy models 
were run for these two exemplars varying a number of parame- 
ters about the adopted solution in Tab. [8] O abundance varied by 
+0.2 dex; stellar black body temperature by 10-15% (20000K 
for F56#2, with T BB 155 000K and 10000 K for F34#l, with 
Tbb 71 000 K); the density by a factor of two; and a model at- 
mosphere rather than a black body. For each modified parameter 
set, Cloudy models were run to match the spectra within errors as 
far as possible by freely varying all the other, non-displaced, pa- 
rameters. Even this limited analysis involved running more than 
180 separate Cloudy models. 

The results of the comparison of Cloudy models are sum- 
marised in Tables [9] and [TO] The third columns list the dered- 
dened spectrum from Tab. Q~2] and the 4th column the spectrum 
resulting from the adopted model given in Tab. [8] The successive 
columns list the resulting spectra from the incremented param- 
eters. The lower half of each table lists the stellar temperature, 
density and abundances corresponding to the spectrum in the up- 
per part of the same column; values in column 4 being identical 
to those in Tab. [8] An overall merit factor, FoM, of the match of 
the model spectrum with the observed spectrum, taking account 
of the S/N of the fluxes, is defined by: 

n 

FoM = J] \{{F obs - F mod )IF obs )\ x (F obs /Err obs )) (1) 
(=i 

where F obs and Err obs are respectively the dereddened flux and 
error and F mo d the model flux; the summation is taken over the 
number of emission lines. Since some lines are in fixed ratio, 
such as H lines, a subset of the lines was adopted to dispense 
with some redundant information (except for the more important 



lines such as Ha and H/3, the weak He I lines and the [O III] 
4959A line). 

Table [9] presents the comparative models for the bright PN 
F56#2 (5601 in Hui et al. 1993a). The comparison of the ob- 
served and model spectra makes it clear that the He I 5876A line 
is strongly under-estimated, most probably because of the poor 
removal of the strong Na I telluric emission. In general satis- 
factory models could be found, except for the cases of O - 0.2 
abundance (column 6 of Tabled and Tbb - 20000 K (column 
8). What constituted a satisfactory match to the spectrum was 
not simply given by the FoM value, but certain line ratios, which 
act as important diagnostics of nebular conditions, e.g. He ll/Hfi, 
[O III] 5007A/H/3 and [O III] 5007/4363A), were given higher 
weight in assessing the quality of the match. It is clear that a 
range of O +0.2 produces model spectra in which the important 
[O III] lines are not well-matched, and thus representative er- 
rors on O in the range 0.1 to 0.15 are suggested. Similarly the 
range of stellar temperature of +20 000 K and density varying 
by a factor 2 also over-estimate the allowed range of some of the 
line ratios. However it is clear from these parameter ranges that 
the error bars are not necessarily symmetric. H, He and PG1 159 
atmospheres from Rauch et al. 12003 1 were input to Cloudy and 
although the choice of temperatures is not very extensive, a sat- 
isfactory fit with the PG1159 atmosphere for 170000 K was 
found. However this should not be considered as a best estimate 
since a comprehensive range of temperature and gravity was not 
tested. 

The comparative models for the fainter PN F34#l (2906 in 
Hui et al. 1993a] W5oo7a=24.96) are given in Tab. [TOl Here the 
number of lines is less than F56#2, since [O III]4363A and no 
He I lines were detected and He II is a limit. The [O II] 7325A 
line is clearly under-estimated in the measured spectrum by com- 
parison with the model spectra. Satisfactory matches could gen- 
erally be obtained with the ranges of parameters listed. The mod- 
els with a range of O of ±0.2 is clearly accommodated by the 
data. Only the lowering of the stellar temperature by 10000 
K (column 8 of fTQb produced an unsatisfactory fit. A PG1159 
model atmosphere (Rauch 2003) of 70 000K produced a good fit 
(column 11), but the 70000 K He model produced a similar fit 
with an O abundance lower by 0.40 dex than the black body. 

From this comparison of observed and model spectra with 
a higher and a lower S/N, we conclude that the typical depen- 
dence of errors on the choice of Cloudy models are generally in 
the range of 0.15 - 0.20 for the O abundance, around 10-15% 
for He and around 0.2 for Ne and N. These values are similar to 
those found by Magrini et al. 2004 comparing abundances from 
Cloudy models with those of the ICF method. These represen- 
tative abundance errors are shown in plots, such as Figs. |6]and 

m 



6. Discussion 

6. 1 . Individual PN spectra and abundances 

A sample of 51 PNe in NGC 5128 have been observed over 
a range of projected galactocentric distances from 1.7 to 18.9 
kpc and covering the PN luminosity function from the bright- 
est known PN in NGC 5128 (PN 5601, m 5 007A = 23.51) to 4.1 
mags fainter. The mean 5007A/H/3 ratio is 11.00 ± 0.55 (r.m.s 
on the mean) for the 42 spectra with the highest S/N. The ex- 
trema of the 5007A/H/3 values for this subset are 6.0 (+ 0.67) to 
18.8 (+3.7); the value of this ratio is very sensitive to the sub- 
traction of the underlying continuum and the variation of the H/3 
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Table 9. Comparative Cloudy photoionization modelling of F56#2 



Species 


A (A) 


Flux& 


Model 


Model 


Model 


Model 


Model 


Model 


Model 


Model 






error 




O + 0.2 


O-0.2 


T sfl + 20 kK 


T ss - 20 kK 


N, x 2.0 


N e x 0.5 


Mod.atmos. 


[Oil] 


3727 


56 ± 13 


59 


58 


47 


53 


59 


32 


53 


57 


[Ne III] 


3868 


124 ± 13 


124 


127 


128 


122 


123 


126 


125 


125 


[0 III] 


4363 


28 ± 8 


24 


15 


27 


26 


22 


22 


22 


27 


He I 


4471 


6± 3 


6 


6 


6 


3 


11 


6 


9 


6 


HeH 


4686 


39 ± 8 


39 


38 


34 


46 


30 


37 


30 


38 


HjS 


4861 


100 ± 


100 


100 


100 


100 


100 


100 


100 


100 


[O III] 


4959 


457 ±31 


449 


448 


398 


451 


471 


454 


458 


456 


[O III] 


5007 


1367 + 90 


1350 


1348 


1199 


1358 


1417 


1367 


1377 


1371 


He I 


5876 


6± 3 


18 


18 


19 


11 


33 


19 


28 


20 


Ha 


6562 


281 ± 19 


286 


285 


290 


286 


286 


284 


285 


288 


[Nil] 


6583 


55+5 


54 


57 


58 


53 


53 


54 


53 


57 


He I 


6678 


7 ± 4 


4 


5 


4 


2 


8 


4 


7 


5 


[Ar III] 


7133 


8± 3 


9 


6 


10 


8 


9 


7 


9 


8 


[O II] 


7325 


14+4 


16 


14 


13 


19 


14 


15 


9 


16 


Model input 






















T Bfi (kK) 






155 


152.5 


150 


175 


135 


150 


132.5 


170 1 


N H (cm- 3 ) 






10000 


10000 


10000 


12000 


8000 


20000 


5000 


10000 


12+log(He) 






11.15 


11.15 


11.20 


11.00 


11.35 


11.15 


11.30 


11.18 


12+log(C) 






8.40 


9.20 


6.60 


8.40 


8.40 


8.70 


8.40 


8.20 


12+log(N) 






7.85 


8.02 


7.84 


7.85 


7.88 


8.01 


8.05 


7.85 


12+log(0) 






8.43 


8.63 


8.23 


8.48 


8.43 


8.49 


8.40 


8.39 


12+log(Ne) 






7.72 


7.96 


7.57 


7.75 


7.70 


7.78 


7.70 


7.67 


12+log(Ar) 






6.00 


6.00 


6.15 


6.00 


6.00 


6.00 


6.00 


5.95 


FoM 






7.53 


8.30 


13.46 


7.34 


14.79 


9.06 


12.90 


8.56 



Notes. 1 Model atmosphere: PG1 159 atmosphere with T« = 170 000 K and log g = 7.00 from Rauch 2003 
S abundance fixed at 12+log(S/H) = 6.55 



absorption line along the slitlet. The mean logarithmic extinc- 
tion correction at H/3 is 0.56 ± 0.20. The overall Galactic ex- 
tinction to NGC 5128 is E fl ^ v =0.123 (Burstein & Heiles [l984l i 
or E B „ y =0.115 (Schlegel et al l 19981 1, the latter is equiv alent to 
c=0.17 for a Galactic extinction law with R=3.2 (Seaton [1979l) . 
The lowest values of extinction measured in the PNe (Tab. [T2l 
are compatible with no local extinction in NGC 5128. The his- 
togram of c values peaks at 0.45 (Eb-v - 0.30) and shows no 
obvious trend with radial offset beyond 200" from the nucleus; 
the five highest values occur within a radius of 200" (3.7 kpc). 
The frequency of higher c values at lower radii is probably due 
to the high line of sight extinction in the vicinity of the dust 
lane, rather than any intrinsic PN property (e.g. high intrinsic 
dust content associated with young PNe, such as NGC 7027; 
Zhang et al. 120031) . 

The ratio of 5007 A/H/? does not show any significant over- 
all gradient with projected radial offset from the galaxy centre; 
the highest values lie in field F42 where the stellar continuum 
is strongest and removal of the underlying H/3 is the most crit- 
ical. There is no evidence for a gradient in [Ne III]/[0 III] im- 
plying a constant ++ /Ne ++ ratio. If there is any O enrichment, 
e.g. during the third dredge-up (Pequignot et al. 2000 1, or deple- 
tion (through the ON cycle, Leisy & Dennefeld 2006) occurring 
in the nebular envelopes of the PNe, then this constant value 
of the ionic ratio implies that any alteration in O abundance 
must be accompanied by a corresponding change in Ne. Since 
O-Ne correlation is not predicted by AGB evolution (Karakas & 
Lattanzio 2003), the rather constant O/Ne ratio implies that O 
enrichment/depletion is not an important effect. 

The mean 12 + log(0/H) is 8.52 + 0.03 (median value 8.48) 
and the range of [O/H] is -0.66 to +0.31. The mean O/H abun- 



dances inside and outside 10 kpc are identical. Of the three PN 
with the lowest O/H, two occur close to the nucleus (F42#4 and 
F42#10) and the other is at 19.5 kpc (F34#2) and all have low 
[O III]/Hj0, without notably strong [O II]; F34#2 has lower N/H. 
The PN with the highest O/H occur in fields 42 and 56 and are 
distinguished by very high [O III]5007/Hy6 but with large er- 
rors. The problems in removing the underlying stellar contin- 
uum from slit measurements can contribute to large uncertainty 
on the line fluxes for fainter nebulae in the central regions (see 
the discussion in Roth (120061 1). and very high [O III]/Hy8 of > 15 
must be viewed with caution since the H/3 may be considerably 
underestimated. 

The mean N/O ratio is 0.51 ± 0.06. Four PNe were found 
with high N/O ratio (F56#8, F56#12b, F56#18 and F42#4) but 
only F56#12b and F56#18 have He/H and N/O which c ould b e 
classified as characteristic of Type I nebulae (Peimbert 119781 1. 
Since F56#12b was modelled by a lower temperature central star 
(the lowest of all the PN modelled, Tab. [8]), it may be that opti- 
cal depth effects and the assumption of a black body may yield 
misleading results. F56#18 is perhaps a better candidate for a 
Type I nebula, having a relatively high stellar temperature and a 
higher reddening, but this PN does not possess the highest lumi- 
nosity central star. The brightest PN, and incidentally the nebula 
with the highest luminosity star (Tab. [SJ, has N/O of 0.12 and 
shows no elevation of this ratio, compared for example to the 
mean Milky Way value of 0.47 (Kingsburgh & Barlow 1994). 
The shorter timescales for PN evolution of Type I PN make them 
minor contributors to a PN population luminosity function, but 
they contribute about ~10% of Milky Way PN by number. 

The mean oxygen abundances in NGC 5128 PNe are inter- 
mediate between values for the LMC (12 + log(0/H) mean 8.4 
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Table 10. Comparative Cloudy photoionization modelling of F34#l 



Species 


^(A) 


Flux& 


Model 


Model 


Model 


Model 


Model 


Model 


Model 


Model 






error 




O + 0.2 


O-0.2 


Tub + 10 kK 


Tbb - 10 kK 


N„ x 2.0 


N e x 0.5 


Mod.atmos. 


[Oil] 


3727 


86 ± 9 


89 


85 


88 


87 


60 


49 


144 


85 


[Ne III] 


3868 


93+6 


91 


92 


93 


92 


95 


92 


94 


96 


Hell 


4686 


2+ 2 


2 


2 


4 


4 


1 


2 


4 





HyS 


4861 


100 ± 


100 


100 


100 


100 


100 


100 


100 


100 


[0 III] 


4959 


393 + 15 


384 


391 


384 


393 


248 


385 


387 


390 


[O III] 


5007 


1164 + 41 


1156 


1177 


1157 


1183 


747 


1157 


1166 


1174 


Ha 


6562 


283 + 12 


288 


290 


287 


286 


286 


287 


288 


287 


[Nil] 


6583 


104 + 6 


104 


104 


104 


104 


103 


104 


105 


103 


[Ar III] 


7133 


4± 3 


3 


4 


4 


4 


3 


3 


4 


4 


[Oil] 


7325 


4± 3 


17 


16 


17 


17 


8 


16 


17 


14 


Model input 






















T ss (kK) 






71 


68 


78 


81 


61 


67.5 


80 


70 1 


N H (cm- 3 ) 






8500 


10000 


7500 


8500 


8500 


17000 


4250 


7000 


12+log(He/H) 






11.06 


11.06 


11.06 


10.90 


11.30 


11.06 


10.90 


11.20 


12+log(C) 






8.60 


8.50 


8.40 


8.80 


8.50 


8.60 


8.30 


8.20 


12+log(N) 






8.25 


8.38 


8.13 


8.16 


8.50 


8.40 


8.03 


8.18 


12+log(0) 






8.73 


8.93 


8.53 


8.60 


9.00 


8.78 


8.54 


8.60 


12+log(Ne) 






8.00 


8.23 


7.78 


7.84 


8.60 


8.07 


7.77 


7.88 


12+log(Ar) 






5.65 


5.85 


5.55 


5.60 


5.90 


6.65 


5.55 


5.60 


FoM 






6.44 


5.79 


6.48 


6.74 


26.62 


10.20 


12.90 


5.30 



Notes. 1 Model atmosphere: PG1 159 atmosphere with T» = 70 000 K and log g = 7.00 from Rauch 2003 
S abundance fixed at 12+log(S/H) = 6.55 



(e.g. Leisy & Dennefeld2006), and that of the Milky Way (mean 
8.68, Kingsburgh & Barlow 1994). A plot of the dependence of 
the [O III]5007A luminosity on 12+log(0/H) (Fig. [5} shows a 
tendency to increase with the oxygen abundance, in line with the 
theoretical relation of of Dopita et al. (1992 ) as observationally 
transformed by Richer (1993 , see Figure 2 and equation 1). A 
second order fit to the observed points demonstrates a compara- 
ble dependence of [O III] luminosity on (O/H) to the theoretical 
relation and is shown in Fig.[5]by a dashed line. 

The region summed spectra provide good detections of a sin- 
gle ionization species of S and Ar (Tab. [6]) and lines of these 
species were detected in a number of the brightest PNe (Tab. 
[TT1and[T2l allowing useful indications of S/H and Ar/H (Tab. [7] 
and Tab. [8]). S and Ar abundances serve an important purpose in 
that they are not considered to be affected by the AGB nucle- 
osynthesis (e.g. Herwig 2005). Whilst O is the element which is 
most easily measured in PN since it is the dominant coolant, it 
can be destroyed by CNO cycling or synthesized during helium 
burning which impacts its usefulness as a metallicity indicator. 
Pequignot et al. (2000) have found evidence for third dredge- 
up O enhancement in a PN in the Sagittarius dwarf galaxy 
(12+log(0/H) = 8.36), but only at the level of <0.03 dex. By 
comparing the lock-step dependence of Ne and O in PN samples, 
Richer & McCall (2008 ) conclude, in agreement with Karakas 
& Lattanzio (2003 ), that the majority of PNe do not show any 
changes in the O or Ne abundances across the AGB and PN 
transition. The Ne/H abundances are well correlated with the O 
abundances for the NGC 5128 PNe as shown by Fig. [6] (slope 
1 .22+0.09, or 1 .17+0.04 excluding the 4 points furthest from the 
linear relation) This result is similar to the value found by Leisy 
& Dennefeld (2006 ) for the LMC (slope 1.13). The relation be- 
tween N/O and O/H is also shown in Fig. [6] There is considerable 
scatter with a weak but insignificant trend for an anti-correlation, 
as also found by Magrini et al. l2004l for PNe in M33. The mean 



36.5 



H 35.5 




8.6 8.8 
12 + log (O/H) 



Fig. 5. The dependence of the oxygen abundance of the 
NGC 5128 PN sample from the Cloudy photoionization mod- 
elling (Tab. [8) on the observed [O III] 5007 A luminosity (uncor- 
rected for reddening) from the 5007 A photometry of Hui et al. 
1993b The dotted line shows the theoretical relation of Dopita 
et al. ( 1 19921 ) scaled to the brightest PN, while the dashed line 
shows a second order polynomial fit to all the data points. 



Ar/H abundance (6.3 for 21 PNe) is indistinguishable from the 
mean for 70 PN in the Milky Way (Kingsburgh & Barlow [1994] ). 
The mean S abundance (6.8 for 9 PNe) also appears to be similar 
to the mean value from the Kingsburgh & Barlow ( 1994) sam- 
ple. One is led to the surprising conclusion that most of the PN 
abundances in NGC 5128 are not significantly different from the 
mean values for Milky Way PNe. 
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Fig. 6. The relationship between log (Ne/H) and log (O/H) is 
shown (upper) and for Log(N/0) against log (O/H) (lower) for 
the individual NGC 5128 PN abundances tabulated in Tab. [8] 
The representative error bars are indicated. 



6.2. PN5601 

PN 5601 (F56#2) is the brightest known PN and the best-studied 
individual PN in NGC 5128, and in deed o f all PN known be- 
yond the Local Group. Walsh et al. (Q999) measured long slit 
observations of this PN amongst a few others and the major dif- 
ference between the spectrum presented in Tab. QT|is the Ha/Hfi 
ratio and the higher He II/H/? ratio. The earlier observations were 
taken with a fixed slit position angle over a considerable range 
of parallactic angle and suffer from wavelength dependent flux 
loss from the slit. This is demonstrated by the very low value of 
extinction derived. The value of extinction derived for this neb- 
ula is now seen to be quite large, but only slightly higher than the 
mean value for all the observed PN (Efi_y=0.32), giving no evi- 
dence that it is intrinsically dustier than the fainter PNe. There is 
some expectation from examples of Milky Way PNe, that young, 
high mass, optically thic k PN ( c.f. NGC 7027, Zhang et al. l2005l 
NGC 6302 Wright et al. 20 11) have higher intrinsic dust extinc- 
tion. 



Fig. 7. The radial dependence of the oxygen abundance (12 + 
log(0/H)) from the observed PN modelled by Cloudy (Tab. [8). 
The points are coded for the three regions as in Fig. F56 by 
open stars (blue); F42 as open squares (red); F34 by filled circles 
(green). A representative error bar on the O abundance is shown. 



6.3. Abundance gradient from PNe 

One of the primary aims of this study was to investigate any 
abundance gradient of the a-elements, as revealed by the PN 
spectra, and to compare it to gradients in the stellar properties 
such as the metallicity distribution function. A plot of the vari- 
ation of O/H v. projected radial offset (Pig. [7j, for the PNe with 
O/H values from the Cloudy models (Tab. [8}, shows two particu- 
lar features. First, the range of O/H values is similar in the inner 
regions (projected radius, % <4 kpc) to the range in the outer 
regions (>13 kpc). Second, there is no evidence of a gradient in 
O/H. A least squares fit confirms this; the linear correlation co- 
efficient is -0.07 for 40 points. Rather, there is a common mean 
value of 8.52 for the inner and outer points (considered as less 
than and greater than % of 10 kpc). However, the highest O/H 
values do occur at % < 4 kpc, in the inner field F42. Although 
our sample is not large, the result that there is no significant O/H 
gradient in NGC 5128 appears to be robust. 

The distribution of [a/Fe] v. [Fe] has an important role in de- 
termining the enrichment history of a galaxy (e.g. Matteucci & 
Recchi l2001l) . O, as representing the a elements, is primarily a 
core collapse SN product released early in star formation activ- 
ity, whilst Fe, as the chief product of Type la SNe, is released 
over many Gyr. Figure [8] shows the Fe v. Mg b EW plot for the 
integrated spectra of regions F42 and F56 with several tracks 
from Thomas et al. (2003) models for selected metallicity and 
[a/Fe] ratio. To be compatible with Harris et al. (2000), who cal- 
ibrated their metallicities by assuming that their stars were as old 
as those in old Galactic globular clusters, then the models must 
be compared at an age of 15 Gyr. In this case the datapoints for 
fields F42 and F56 lie close to the tracks of [a/Fe] = 0.3, giv- 
ing an [a/Fe] ratio of 0.25, and show metallicities [Z] around - 
0.6 (there are no tracks at metallicity intermediate between -0.33 
and -1.35 in Thomas et al. l20031 >. If, however, the age of the stars 
were younger (as indicated e.g. by Rejkuba et al. (2005. 2011)), 
the derived metallicities would be about 0.2 higher, giving an 
average value of -0.4. This would imply an oxygen abundance 
relative to solar of -0. 15, a value that is consistent with the mean 
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O abundance of -0. 17 with respect to the solar value (Scott et al. 
120091 . 

These abundances can also be compared with the MDF from 
Harris & Harris (12002) in three fields at 7.6 kpc SW of the cen- 
tre, at 20 kpc S and 29 kpc S (distances rescaled to 3.8 Mpc). 
Fig. [9] shows the histogram of all the PN O/H abundances with 
respect to Solar overplotted against the completeness corrected 
MDFs in the inner field (7.6 kpc) and two combined outer fields 
(Table 1 of Harris & Harris 2002). The PN sample covers the 
range around the inner and outer fields of Harris et al. (|2002 1 . 
A shift of -0.25 is required to align the peak of the PN O/H dis- 
tribution with the stellar MDF in the inner field (magenta his- 
togram). However it is strongly apparent that the distribution of 
(O/H) for the PNe is much narrower than for the stellar metal- 
licity, so matching the distribution peaks may not be a justifiable 
simplification. Splitting the PN into two samples with positions 
greater and less than 10 kpc shows a shift in the peak value of 
(O/H) of — 0.2 for the inner field and — 0.3 for the outer field, al- 
though the number of PNe with abundances is lower in the outer 
field (11). For the full PN sample, with an average O abundance 
of -0.17 compared to Solar, the total metallicity is then -0.4 as- 
suming a mean [a/Fe]=0.25. Harris et al. (2002) find an average 
metallicity of -0.6 in their outer fields, which becomes -0.4 as- 
suming an age of 8 Gyr; thus taken at face value, assuming the 
metallicity of the field stars is that of the PN, the age of these 
stars in NGC 5 1 28 is around 1 Gyr. 

Both methods - matching observed spectral indices to stellar 
population models and comparison with photometrically derived 
stellar metallicity distributions - provide an average [O/Fe] ratio 
of 0.25 at a stellar metallicity (viz. [Fe/H]) of around -0.4, and 
an average age of 8 Gyr. This single value is obviously a sim- 
plification for the behaviour in a whole galaxy, but demonstrates 
that the [O/Fe] v. [Fe/H] value for NGC 5128 is within the range 
of a-enhanced stellar metallicities for early type galaxies (c.f. 
Thomas et al. 2003). The modelling of the colour-magnitude di- 
agram by Rejkuba et al. (1201 11 1 also independently favours a- 
enhanced stellar abundances in NGC 5128. 

However it can be argued that while all the stars are repre- 
sented in the integrated stellar spectra and photometric surveys, 
only a subset of these stars become PNe. So, the mean values 
of the two samples of the PN [O/H] and the stellar [Fe/H] in 
Fig. [9] are not necessarily strictly comparable. In particular the 
PNe may derive from a younger population, on average. Only 
the higher Z stars produce observable PN, while the stellar con- 
tinuum is comprised of stars of all metallicities. Overall the oc- 
currence of PN is highly delimited relative to all the stars in a 
galaxy and in addition the PN central star masses also have a 
very restricted range relative to all white dwarfs. Thus the PNe 
may only come from a young population; that would explain the 
tendency toward both higher Z (Fig. |9) and higher central star 
masses (see Fig.fTOl. The old stars, which have lower Z on aver- 
age, fail to produce PNe because their masses are too low. 

Chiappini et al. d20091 > present a careful comparison of the 
spectroscopically determined abundances of Milky Way Bulge 
stars with Bulge planetary nebulae. Their sample of 166 PN in 
the Bulge with well-determined abundances shows a broad dis- 
tribution of 12 + log(0/H) from about 8.2 to 8.9, with a mean 
of 8.57, rather similar to the mean value and range determined 
here in NGC 5128. However the O/H abundance collated from 
various spectroscopic studies of Bulge giant stars shows a lower 
value by ~0.3 dex. Various suggestions are offered by Chiappini 
et al. (120091) for this discrepancy, but none appears to be con- 
clusive. Comparing O from the Bulge PN with [O/Fe] from the 
stars, implies [OpN/Festars] ~ -0.1. This value is not strictly com- 
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Fig. 8. The Lick index < Fe > is plotted against Mg b for the 
integrated spectra of the two inner fields F42 (red lozenge) and 
F56 (blue lozenge). Two families of tracks from Thomas et al. 
(2003) for [a7Fe] of 0.3 (open symbols) and 0.0 (filled symbols) 
are shown for metallicity [Z] of -1.35, -0.33, 0.0 and 0.35 by 
lines with points. The lower set of four lines have [a/Fe] = 0.3, 
and the upper ones [a/Fe]=0. Models with [Z]=0.35 are indi- 
cated in black, [Z]=0.0 in orange, [Z]=-0.33 in green and [Z]=- 
1.35 in magenta. The ages in Gyr along each track are indicated. 



parable to that in NGC 5128 since [Fe/H] is derived from pho- 
tometry; spectroscopy of individual stars in NGC 5128 is ruled 
out until larger telescopes are available. Taken at face value the 
difference in [O/Fe] could imply that the formation history of 
NGC 5128 as a classical giant elliptical, is dissimilar to that of 
the Bulge; however the discrepancy in O abundances of Bulge 
PN and giants found by Chiappini et al. (2009) weakens this 
conclusion as does the large spread in stellar Fe abundances (in 
both NGC 5128 and the Bulge). 

6.4. PN and the star formation history 

One of the advantages of running Cloudy models to determine 
the PN abundances is that the central star luminosities and tem- 
peratures are derived. Since the PNe are optically thick, the 
match to the spectrum effectively provides the temperature and 
the m 5007 ^ photometry, combined with the reddening, provides 
the luminosity. Fig.[10]shows the Log L, Log T e f/ points for the 
40 PNe modelled (Tab. [8]). The log L values span a range of a 
factor 6, but are, not surprizingly, at the high end for PNe central 
stars, since only the top few magnitudes of the PNLF is explored 
by these data. No radial trend in the central star luminosities is 
apparent, paralleling the lack of an O abundance gradient (Fig. 
13. A few evolutionary tracks from Vassiliadis & Wood (1997) 
for higher mass progenitor stars are shown on Fig. [10] although 
there is currently no information on whether the central stars are 
on the H or He burning track. For the H burning tracks the bulk 
of the points match stars with masses above 2 M Q . 



16 



Walsh, J. R., et al.: Abundance Distribution in NGC 5128 from PNe 



MDF 20 + 30 kpc 

MDF 8 kpc 

PN 



-1.5 



-0.5 
[O/H] / [Z] 




Log T.„ (K) 



Fig. 9. The histogram of the log O/H abundance relative to Solar 
is plotted for all PN with O/H determined from Cloudy models 
(Tab. [8) as a continuous line (black). For comparison the stel- 
lar metallicity distribution function (MDF) from Harris & Harris 
(2002) is plotted in two regions: an inner region (8 kpc) dotted 
(magenta); a composite of two outer regions at 20 and 30 kpc 
dashed (green). 



There are multiple routes in terms of stellar mass, metallicity 
and mode of turn-off from the AGB (H or He burning track) to 
reach the same PN luminosity, so relating a given PN in a galaxy 
to its progenitor star formation episode is non-trivial. Marigo 
et al. d20041 > have modelled effects of stellar populations on the 
[O III] cut-off of a population of PN, as measured for the PNLF 
They find that the peak [O III] luminosity is emitted by stars 
with initial (main sequence) masses of about 2.5 M with a very 
strong dependence with age. PN from higher mass progenitors, 
although being intrinsically luminous, evolve very rapidly and 
so make a minor contribution to the bright end of the PN pop- 
ulation. In addition, Ciardullo & Jacoby (1 19991 1 showed that the 
intrinsically massive and luminous progenitors are prolific dust 
producers, and therefore self-extinct. Thus the most massive pro- 
genitors will generally appear as rare and faint PNe in a galactic 
sample. Ciardullo et al. (2002) suggest that the progenitors of 
the most luminous PNe in galaxies, if they are not formed from 
stars younger than ~ 1 Gyr, may occur predominantly in binary 
systems. 

In principle the calibration of PN luminosity v. stellar pro- 
genitor mass can be used to identify the age of the starburst 
giving rise to an observed PN population. If indeed the bulk 
of the PN observed in this study arise in progenitors of about 
2.5 M Q , then from the grid of evolutionary tracks of low mass 
stars of Girardi et al. (|2000l one can infer a progenitor age of 
0.7 Gyr (see Fig. [TT] derived from the Girardi et al. tracks for 
[Z/Z ] 0.0, -0.4 and -0.7 and the stellar age to AGB termination). 
Such an age would place the PN in a very recent star formation 
episode. This seems rather unlikely since there is little evidence 
for an extensive stellar population so young. A young (2-4 Gyr) 
minority (20-30%) component was inferred from the HST stel- 
lar photometric studies of Rejkuba et al. (120051 and [2011 i, but 
not with stars younger than 2 Gyr as appears to be required to 
match the PN. Younger stars are brighter and so should have 
been well-detected in optical photometry, ruling out their pres- 
ence in NGC 5128. 



Fig. 10. An HR diagram showing the black body log L v. log T e ff 
coordinates of the central stars of the PNe with spectra modelled 
by Cloudy (Tab.[8]l and a selection of low mass H and He burning 
model tracks from Vassiliadis & Wood (11997) . The points are 
coded (and coloured) for the three regions (F42, F56 and F34) 
as in Fig. [2] Representative error bars on temperature at higher 
and lower temperatures are shown; the error bars in Log L were 
computed from the errors on dereddened absolute H/3 flux in Tab. 
M 
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Fig. 11. The age at which a star terminates the AGB is shown 
as a function of the initial main sequence stellar mass for three 
metallicities of 0.019 (long-dashed line, red), 0.008 (Solar metal- 
licity; continous line, black) and 0.004 (approximately LMC 
metallicity; short-dashed line, blue). The curves are derived from 
the stellar evolution tracks of Girardi et al. (20001. 



If the PNe studied here belong to this minority younger pop- 
ulation of age 2-4 Gyr, then the deduced stellar masses, as indi- 
cated by the HR diagram of the PN central stars in Fig. [10] are 
still unexpectedly high. However since these PN belong to the 
bright end of the luminosity function, there is a selection effect 
in favour of the more luminous (viz. highest mass) PN progenitor 
stars. If however the observed PNe arose from an intermediate 
age population (5-8 Gyr), the progenitor stars would only have 
a mass of ~ 1 .2 M Q and the discrepancy between their deduced 
masses and age is more difficult to reconcile with evolutionary 



17 



Walsh, J. R., et al.: Abundance Distribution in NGC 5128 from PNe 



tracks of single low mass stars (e.g. Vassiliadis & Wood 1997]'. It 
would be profitable to study spectroscopically a group of lower 
luminosity PNe in NGC 5128 to determine if their properties 
differ from the high L PNe presented in this study. 

7. Conclusions 

Low resolution spectroscopy of 51 PN in NGC 5128 at a range 
of galactocentric distance of 2-20 kpc have been obtained with 
the ESO VLT and FORS1 instrument in multi-slit mode (MOS) 
in three fields. The PN are drawn from the upper 4 magnitudes 
of the PN [O III] luminosity function. The emission line spectra 
have been analysed and lines typical of ionization by hot PN cen- 
tral stars have been measured. The weak [O DTJ4363A line was 
just detected in 20% of the PNe. In order to determine element 
abundances for a larger fraction of the observed PNe, photoion- 
ization modelling was conducted with Cloudy for 40 PN with 
the highest quality spectra (representing the upper 2 mag. of the 
PNLF). He, N, O and Ne abundances were determined for all the 
PNe, and S and Ar for about half of the sample. For the most re- 
liably estimated element, oxygen, no radial gradient is seen and 
the slope of O/H v. projected radius is identical inside and out- 
side of 10 kpc. The range of [O/H] in the sample spanned -0.66 
to +0.31 dex with the mean 12 + log O/H of 8.52 (median 8.48). 

The PN O abundances were compared to the stellar abun- 
dance measured from Lick indices from the continuum spec- 
troscopy (i.e. along the MOS slitlets not occupied by PN emis- 
sion) and from resolved star photometry from the literature . If 
the stars in NGC 5128 have an average age of 8 Gyr, the stellar 
and PN metallicities agree in the outer parts of NGC 5128 with 
average values around -0.4 Solar. The deduced masses of the PN 
central stars implies progenitor masses above 2 M , favouring 
their formation from a very young component of the interme- 
diate age stellar population with age t5 Gyr. This discrepancy 
between the age of the stellar population and the mass of the PN 
progenitor stars is in line with other studies of the high mass end 
of the PN luminosity function. 
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Table 11. NGC 5128 PN spectra - Observed line fluxes 
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Table 12. NGC 5128 PN spectra - extinction corrected line fluxes 
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